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ABSTRACT 

We analyze archival iKT/STIS/FUV-MAMA imaging and spectroscopy of 13 com- 
pact star clusters within the circumnuclear starburst region of M83, the closest such 
example. We compare the observed spectra with semi-empirical models, which are 
based on an empirical library of Galactic and B stars observed with IUE, and with 
theoretical models, which are based on a new theoretical UV library of hot massive stars 
computed with WM-Basic. The models were generated with Starburst99 for metallic- 
ities of Z = 0.020 and Z = 0.040, and for stellar IMFs with upper mass limits of 10, 
30, 50, and 100 M . We estimate the ages and masses of the clusters from the best fit 
model spectra, and find that the ages derived from the semi-empirical and theoretical 
models agree within a factor of 1.2 on average. A comparison of the spectroscopic age 
estimates with values derived from .ffi)T/WFC3/UVIS multi-band photometry shows a 
similar level of agreement for all but one cluster. The clusters have a range of ages from 
about 3 to 20Myr, and do not appear to have an age gradient along M83's starburst. 
Clusters with strong P-Cygni profiles have masses of a fewxl0 4 MQ, seem to have 
formed stars more massive than 30 M Q , and are consistent with a Kroupa IMF from 
0.1 — 100 Mq. Field regions in the starburst lack P-Cygni profiles and are dominated 
by B stars. 

Subject headings: galaxies: individual (NGC 5236) — galaxies: nuclei — galaxies: 
starburst — galaxies: star clusters — galaxies: stellar content 



1. INTRODUCTION 



In the context of disk galaxies, a nuclear or circumnuclear ring is a region of l arge gas surface 
density and enhanced star formation, located within 2 kpc of the galactic nucleus. I Comer on et al. 
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(|2010l ) distinguish between nuclear disks and nuclear rings, by setting the maximum width of 
the ring to half the ring radius. The most compact ri ngs (ultra-compact nuclear rings, UCNRs) 
have radii smaller th an 200 pc (jComeron et al.l l2008l ). The Atlas of Images of Nuclear Rings 
(jComeron et al.ll2010l ) is the most complete atlas of galactic nuclear rings to date and includes 101 
barred and unbarred disk galaxies with rings. The dust nuclear rings around elliptical and SO galax- 
ies constitute a different class of galactic ring with no associated star formation. The present work 
i s not concerned with th e latter rings. Circumnuclear rings in disk galaxies give rise to starbursts 
(jSersic fc Pastorizalll965l ; lMaoz et al.lll996l ; iBonatto et al.lll998l ). which observati onally are ring- or 
arc-shaped. Note that the term "starburst" does not have a general definition (jConti et al.l 12008k 
Knapen fc Jamesl l2009i ) . A "circumnuclear starburst" is a region composed of star condensations 
with individual masses ranging from a few xlO 4 to greater than 10 6 M Q , bright in the ultraviolet 
(UV) because they have recently emerged from their birth clouds (the nuclear galactic ring), and 
contain hot and massive O and/or B stars. 

Proposed scenarios for the origin of star-forming galactic rings are galaxy collisions or merg- 
ers, accretion of intergalactic gas, and resonances caused by the ac tions of rotating bars o r other 
non-axisymmetric disturbances on the motions of disk gas clouds (jButa h Combed 119961). Note 
that i n their analysis of Ha and continuum images of 327 local disk galaxies, iKnapen fe James 
(|2009l ) found no significant increase in the central concentration of star formation as a result of 
the presence of a close companion. On the other hand, th ere is strong observat ional evidence 
for a connection between bars and circumnuclear starbursts (jHunt fc Malkan 1999). In addition , 
non-starburst barred galaxies may eventually become starburst barred galaxies (jJogee et al.ll2005l ). 
The location of nuclear rings in barred galaxies may be set by inner Lindblad resonances (ILRs), 
which can prevent gas from flowing further i n. Such resonances ca n cause gas to accumulate in 
circumnuclear rings where clusters can form (jButa &: Combed 1 19961). Al t ernat ive m odels for the 
locati on of st arbursts in barred gal a xies c an be found in Shlosman et al. (jl99ol ) and Regan et al 



(j 19971 ) , while Ivan de Ven &: ChaneJ (|2009i ) have a model for the migration of circumnuclear star 
clusters and nuclear rings. 

Circumnucle ar starbursts in b arred galaxies are perhaps the most numerous class of nearby 
starburst regions (|Maoz et alJboOll ). Their detailed study can provide a stepping stone for studies 
of starbursts at higher redshifts. Indeed, lo cal starburst galaxi es (including circumnuclear star 



bursts), show star formation rate densities (Meurer et 



al 



Adelberger fc Steidelll2000i ). and spectral morphologies (IConti et al 



199?]), UV colors JMeurer 



ilar to those of high redshift Lyman-break galaxies (also see 



1996 



Pettini et al 



Heckman et al 



20051 ). 



et al. 1999 



2000), 



sim- 



Determining the ages, masses, and stellar initial mass functions (IMFs) of individual star 
clusters within circumnuclear starbursts is important for understanding the relations between (1) 
galaxy dynamics and interactions and the origin of cir cumnuclear star clusters; (2) ma ssive star 



feedback and the fueling of active galactic nuclei ( AGN. IChen et al 



2009 



Hao et al.ll2009l ): and (3) 



young massive star clusters and globular clusters ( de Grijs 2009 : Elmegreen 201ol ). In this study. 



we determine the ages and masses, and constrain the IMFs of 13 star clusters in Messier 83 (M83, 
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NGC 5236), which hosts the nearest example of a circumnuclear starburst in a barred galaxy. 



M83 is a nea rby f4.5±0.3 Mpc. iThim et al.ll2003r ). southern (Dec[J2000] = -29°51'56"), nearly 
face- on (i = 24°, iRogstad et al.lll974h . SAB(s)c grand-design spiral ga laxy (Ide Vaucouleurs et al 



19921 ). with an optical disk spanning 12'. 9 x 11.5' (17.0 kpc x 15.2 kpc). IWiklind et al.1 (120041 ) have 
estimated that at the distance of M83, a central black hole of moderate size and luminosity would 
be detectable. However, the galaxy shows no indicati on from radio or X-ray observations of hosting 



an AGN (Cowan fc Branch 



the dwarf irregula r galaxy NGC 5253 ( 
40 - 50 pc in size (jCalzetti et al.lll997l . 



19851 ; ISoria Wul 120021) . M8 3 has a close dynamical companion in 



Rogstad et al. ,119741 ) which contains a st arburst nucleus of 



Calzetti et alJll999 . iTremonti et alj|200ll ) 



M83's arc-shaped cir cumnuclear starbur st has been observed at wavelengths ranging from the 
radio to the X-rays (see iDopita et al.1 120101 for references). It spans about 200 pc in length and 
35 pc in thickness, and it is located between tw o rings of dusty gaseou s material which may be 
coincident with two inner Lindblad resonances f El megreen et al.l Il998l ) . The starburst qualifies 
as an UCNR (ring radius of < 200 pc), and as shown in this work, it breaks up into about 20 
compact FUV bright star clusters. Two proposed scenarios for its origi n are, the interaction of 
M83 with its neighbor NGC 5253 about 1 Gyr ago (jvan den Ber gh 1980), and/or the merging of 
M83 with a dwarf satellite in the past. The former scenario is su pported by the "fing er" of H 
I gas projecting outward from the halo of M83 toward NGC 5253 (jRogstad et al.lll974h . and by 
the presence of off-disk H I clou ds in M83 ([Miller et al.ll2009l ). T he merger scenario i s based on 



kinem atic studies using R band (|Diaz et al 



2006b! ) . and submillimeter ((Sakamoto et al 



2006a), near- infrared (jThatte et al. 



2000 



Diaz et al 



20041 ) spectroscopy, which have revealed the presence 



of a controversial second nucleus, hidden by dust, more massive than the optical nucle u s, an d 
located 4" ± 1" northwest of the latter. Recent N-body simulations by iRodrigues et al.l (|2009f ). 
which account for this second nucleus, predict the d i srupt ion of the circumnuclear starburst in less 
than an orbital time. However iHoughton Sz Thattd (|2008l ) have concluded that there is no second 
nucleus in M83, since they cannot confirm th e local maxim u m in the stellar velocity dispersion at 
the location of the "second nucleus" found by iThatte et al.l (|2000l ) . A third scenario for the origin 
of the starburst could be internal galaxy dynamics associated with the bar. 

The ages of circumnuclear star clusters in M83 have been preyiousl y estimated. Based on the 
analysis of CO absorption and Br7 emission- line data lPuxlev et al.l (|1997l ) found that t he northwest 



ern en d of M83's starburst is the youngest, while the southeastern end is the oldest. lHarris et al 



(|200ll ) compared HST /WFPC2 photometry with theoretical population synthesis models and found 
that most of M83's circumnuclear star clusters are less than 10 Myr old, with clusters younger than 
5 Myr preferen t ially l oca ted along the north western end. The statistical analysis of the photometry 
of lHarris et al.l (120011) bvlDfaz et al.1 (l2006bl ) confirms the age gradient along the starburst. In their 
analysis. iDfaz et al.l (|2006bl ) divided the 100° starburst region into 25° angular sectors, and consid- 
ered the mean cluster age as well as the oldest age within each sector, finding 5 Myr for the oldest 
cluster in the northwest and 25 Myr for the oldest cluster in the southeast. The above age gradient 
suggests that star formation occurs in an ordered manner, and supports the idea that the starburst 
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is fed by the inflow of bar-driven material. If real, this gradient has profound implications for the 
formation and evolution of circumnuclear starbursts. For this reason, in this work, we age-date 
13 clusters within the northern portion of M83's starburst, using three different techniques, i.e., 
by comparing archival HST /STIS FUV spectroscopy of the clusters with semi-empirical and new 
fully theoretical model spectra, and by deriving the ages from high spatial resolution HST /WFC3 
multi-band optical photometry. Although our study does not cover the full extent of the starburst's 
arc, this does not affect our results. In addition, we analyze the massive star population in regions 
between the clusters. 

Another goal of this work is to constrain the upper IMFs of our clusters. This is particularly 
interesting in the metal rich circumnuclear environment of M83. M83's circ umnuclear H II regions 



have an oxygen to hydrogen number density ratio of 12+log(O/H)=9.1±0.2 (jBresolin Kennicutt 



2002), and a corresponding fraction of all metals to all elements by mass (hereafter, metallicity) 
that is int ermediate between Z =0.020 and Z=0.040, and supersolar (12+log[O/H] =8.69±0.05, 
Z =O.O13, Asplund et al .112009 ). We estimate the cluster masses from spectroscopy and photometry, 



and check which clusters have formed stars more massive than > 30 M . The random samp ling of 



the stellar IMF may prevent this from happening in low mass clusters (jWeidner et al.ll2010l ). 



This paper is organized as follows: § [2] describes the observations and data reduction, § [3] 
presents the models used in the estimate of the spectroscopic ages and masses of the clusters, § [4] 
explains our procedure for determining the spectroscopic and photometric properties of the clusters, 
in § [5] we present and discuss our results, finally, we conclude in § El 



OBSERVATIONS AND DATA REDUCTION 



We analyzed archival FUV imaging and spectroscopy of compact star clusters and field regions 
located between about 5" (=110 pc) and 9" (=198 pc) of M83's optical center, and spanning an 
arc of ~200 pc in size. The data were taken with the STIS onboard the HST as part of programs 
GTO-7577 and GTO-8785 (PI: Sara Heap). We complemented these data with optical images of 
the same clust e rs tak en with HST's WFC3 as part of program GO/DD-11360 (PI: O'Connell). 
Chandar et al.l (|2010l ) discuss properties derived from optical photometry of star clusters within 
the larger field of view covered by the WFC3 images. Our main focus is on the FUV data, but we 
also briefly discuss the optical data. 

Figure [T] shows a false-color logarithmic-scale image of M83's arc-shaped circumnuclear star- 
burst. The image is a composite of the STIS image taken with the F25QTZ (FUV) filter, and 
WFC3 images taken with the F657N (Ha) and F555W (V) filters. North is at the top and east is 
to the left. The optical nucleus of M83 is the large yellow-green dot in the upper left. The starburst 
starts northwest of the nucleus, wraps about the nucleus towards the southeast, and is punctuated 
by about 20 UV bright clusters. The FUV bright clust ers appear white, w hile fainter ones appear 
blue. The location of SN1968L, which was identified by lDopita et al.1 (|2010l ). is indicated. At M83's 
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distance, 1" represents about 22 pc. The two dusty rings discovered by lElmegreen et al.l ((1998) 
and embracing the starburst are clearly seen. The clusters with available STIS FUV spectroscopy 
are labeled. Note that although clusters 3, 4 and 7 seem to be composed of several bright knots 
clustered together, one knot always dominates the brightness. Figure [1] was kindly produced for us 
by W. P. Blair. 



2.1. FUV Imaging and Spectroscopy with STIS 

The FUV data set consists of a direct image and four two-dimensional (2-D) spectrograms. 
The direct image is projected onto the MAMA detector, which has a field of view of 25"x25" and 
a plate scale of .024 pixel -1 . The spectrograms, one long slit (25"x2") and three multi-object 
(slitless), were taken with the G140L grating, which has an average dispersion of 0.584 A pixel -1 . 
Table [1] gives the observation date (col. 2), exposure time (col. 3), aperture position angle in 
degrees East of North (col. 4), filter name (col. 5), wavelength coverage for the detector + filter 
combination (col. 6), and extension of the most processed data file available from the Hubble Data 
Archive of the FUV observations. 

The data were partially processed through the On-the-Fly Reprocessing (OTFR) pipeline of 
the archive, at retrieval time. This included data linearization, dark subtraction, flat fielding, and 
summing of the repeated observations. In addition, the pipeline corrected the direct image for 
geometric distortion. After processing through the pipeline, in order to increase the signal-to-noise 
(S/N), we added together exposures taken with the same instrument configuration, i.e., the two 
slitless exposures taken with the F25SFR2 aperture, which we refer to as the F25SFR2 exposure. 

Because the spectrograms were taken in either a slitless or a wide slit mode with the G140L 
grating, the spectral traces have large offsets in the dispersion direction. Therefore, the OTFR 
pipeline did not wavelength calibrate or convert these spectra to absolute flux. In order to extract 
wavelength-calibrated one-dimensional (1-D) spectra from the 2-D spectra, we first located the 
position of each spectral trace along the cross dispersion direction, and then used the X1D module 
(McGr ath, Busko, Jk Hod ge, STIS ISR 99-03) available in the Image Reduction and Analysis Facility 



(IRAF. lTodylll986l . Il993l ) for this task. Several iterations were required in order to find the optimal 
value of xoffset, which is the parameter used by X1D for shifting a spectrum along the dispersion 
direction. We used the absorption feature C II AA1334, 1335 + A1336 (hereafter C II 1335), for 
finding the optimal value of xoffset. In our spectra, this feature is unresolved and is a blend of 
Milky Way interstellar and M83 stellar components. 

For the width of the extraction box we used the value recommended in the STIS manual for 



point sources, i.e., 11 pixels (jBrownl 120021 ) . Background spectra extracted from outer portions of 
the 2-D spectrograms were assumed to represent purely geocoronal emission (air glow) from H I 
Lyman-a (Lya) at 1216 A and/or, [O I] 1302 + 1305 + 1306 A, depending on the aperture, and were 
subtracted from the cluster spectra. We did not subtract the underlying diffuse light because of its 
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large variations as a function of position. 



The direct image and the 2-D spectrograms from the STIS are shown in the four panels of 
Fig. [2l Panel (a) features the 25"x25" direct image in logarithmic scale, with the location of M83's 
optical center marked with a cross, and the spatial scale marked with a vector. Panels (b) and 
(c) show the long-slit and the multi-object spectrograms respectively, both taken with the same 
aperture (25MAMA). The vertical band in panel (b) is due to Lya airglow filling the 2" -wide slit (it 
is seen more clearly in the electronic version). The Lya airglow contaminates much of the slitless 
spectrogram shown in panel (c). Finally, panel (d) features the sum of the two slitless spectrograms 
taken with the F25SFR2 aperture, which rejects the Lya airglow. No contamination from the [O I] 
air glow is apparent in panel (d) because the emission brightness geocoronal [O I] is typically less 
than 10% of that of geocoro nal Lya, depending on the t ime of observation and the position of the 
target relative to the Earth (jMiskey fc Bruhweilerll2003l ). 



We used the slitless exposures for extracting spectra of those clusters falling outside the slit, 
as well as for studying the effect of rejecting the Lya air glow on our results. 

The spectral resolution, depends on the spatial extent of the star clusters. In order to estimate 
the spectral resolution of our spectra, we determined the full width at half maximum (FWHM) 
of the cross dispersion profiles of the spectral traces. For our point-like representative, GD 153, 
the FWHM is 0.09". Due to blending with the profiles of close neighbours, we were only able 
to determine the FWHM of four clusters by our method. For clusters 1, 4, 10, and 11+12 we 
obtained FWHM of 0.17", 0.30", 0.19", and 0.23", respectively, which suggests that the clusters 
are partially resolved. Since these are neither point sources nor fully extended sources, we expected 
the spectral resolution of our long-slit spectra to be intermediate between the values given in the 
STIS instrument handbook, i.e., 0.9A (point source) and 48A (extended source). Multi-component 
Gaussian fitting of the C II feature at AA1334, 1335 + A1336, carried out with the Peak Analysis 



code (available at ftp:/ /ftp. ncnr.nist.gov/pub/staff/dimeo/pan. zip) reveals that the FWHM of this 
feature ranges from 2.5-7 A. Our field spectra are expected to have the resolution of a fully extended 



source. 



We corrected the cluster spectra for reddening due to the presence of dust in t he Galaxy and 
in M 83. For the Galaxy's contribution to the reddenin g, we adopt E(B — V) = 0.06 (jSchlegel et al 



19981 ). and we assume the Galactic extinction curve of iFitzpatrickl (|1999l ). For M83's contribution, 
we (i) fitted a power law of the form F oc \^ to the continuum in the range 1240 — 1600 A, (ii) 
assumed that any deviation of j3 from —2.6, which is the expected value for a dust free starburst, is 
due to extinction by local dust, and (in ) dereddened the obs erved spectrum until —2.6 was reached, 
using the starburst obscuration law of lCalzetti et al.1 (120001 ) . 
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2.2. Matching the Clusters with their Spectral Traces 

In order to match each cluster with its FUV spectral trace, we started by directly comparing 
the positions of the clusters in the STIS direct image, with the positions of the spectral traces in the 
2-D spectrograms. For this, we rotated the image and vertically shifted the spectrograms until the 
orientations matched, as shown in Fig. [3l which shows enlarged portions of the STIS image, side 
by side with two spectrograms. Note that the slitless spectrograms all have the same orientation, 
while the orientation of the long-slit spectrogram is different. Therefore, in Fig. [3j we show only 
the slitless spectrogram taken with the F25SFR2 aperture. In panel (a) of Fig. [3j we indicate the 
position of the long-slit, the locations where we extracted the sky and field region spectra, and the 
correspondence between clusters within the slit and their spectral traces. Panel (c), which shows 
the same image as panel (a), but with a slightly different orientation, shows only one sky region 
due to the fact that the second falls outside of the field displayed in the panel. 

Unfortunately, and as can be seen in Fig. [3j in spite of the different orientations of the spec- 
trograms, the spectral traces of the two brightest clusters, 11 and 12, always overlap. Furthermore, 
as seen in panel (d), the faint spectral traces of clusters 4 and 5 also overlap, as do the traces of 
cluster 14 and a source not in our sample. Fortunately, we can use the long-slit spectra for clusters 
4 and 14. However, cluster 5 was rejected from our analysis due to its poor spectral trace. 

The spatial correspondence between star clusters and their spectral traces can be checked in 
Fig. HI where we compare the cluster radial profiles derived from the STIS direct image (bottom 
curves in left and right panels), with the cross dispersion profiles derived from the long-slit (top 
curve, left panel) and the slitless (top curve, right panel) 25MAMA spectrograms. Also shown in 
Fig. [J] is the cross-dispersion profile of a point-like source, represented by white dwarf GD 153 (top 
dotted curve in left and right panels), observed as part of program 7097, in the same configuration 
as the long-slit spectrogram. For convenience, we divided each curve in Fig. H]by the amplitude of 
the tallest peak, which for the clusters is the amplitude of target 11+12 (which overlap), and for 
the white dwarf is its own amplitude. For clarity, the top curves in Fig. [3] (cross dispersion profiles) 
were all shifted by 0.2. The clusters and regions used for the background subtraction are labeled. 

For the radial profiles, we first rotated the image to match the orientation of each spectrogram, 
and we then obtained the count rate at a given vertical position by integrating over the horizontal 
direction. The radial profiles in the left and right panels look slightly different because of the 
different angles by which we rotated the image. We note that although the peaks for clusters in our 
sample are dominated by emission from the clusters themselves, diffuse light and fainter sources 
also contribute to their shape and amplitude. This however, does not affect the identification of 
the spectrum of each cluster. 

For the spectral profiles, the count rate at a given position in the cross dispersion direction 
was obtained by integrating over the dispersion axis. Note the higher background of the top right 
curve with respect to the top left curve, which is due to contamination with Lya air glow over a 
larger portion of the detector in the slitless spectrogram (right panel of Fig. 
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2.3. Optical Imaging with WFC3 

Circular aperture photometry was performed on the clusters using the IRAF task PHOT, using 
an aperture radius of 3 pixels and a background annulus with radii of 10 and 13 pixels. Aperture 
corrections, based on the estimated size of each cluster, were applied to convert the fixed aperture 
magnitudes to total magnitudes. The instrumental magnitudes in the F336W, F438W, F555W, 
F658N, and F814W filters were converted to the VEGAMAG photometric system by applying 



zeropoints available at the following URL: http://www.stsci.edu/hst/wfc3/phot_zp_lbn, We refer 



to the broad-band filters as the "f7" , "B" , il V" , and "J" magnitudes for simplicity, although they 



were not transformed to the Johnson-Cousins system. For more details, refer to IChandar et al 



(|20im . 



3. MODEL SPECTRA 



We compared our observed redden ing corrected FUV spec tra with synthetic dust free spectra 
generated wi th the Windows vers i on (ILeitherer &: C hen 2009) of the widely used package Star- 



burst99 (S99. ILeitherer et al 



1999 



Vazquez Leithererll2005l ). The main input parameters in S99 



are 1) the star-formation law, 2) the IMF, 3) the metallicity and stellar evolution tracks, and 4) 
the stellar spectral library. 



We tried the following parameters. 1) A single burst of star formation, since lBresolin &: Kennicutt 



(|2002l ) ruled out the possibility that star formation within individual clusters proceeded continu- 
ously. Indeed, in M83's circumnuclear region, the He I 5876 absorption line is detected in H II 
regions ionized by the youngest clusters, but not in those ionized by clusters older than about 
6 — 7Myr. If star formation was continuous, the line strength ratio of He I 5876 to H/3 would 
be the equilibrium value in the older clusters. 2) IMFs with hig h mass limits of 10, 30, 50, and 
100 Mq. 3) The stel lar evolution t r acks for non-rotating stars of ISchaller et al.l (|1992l ) for masses 
M < 12 M , and of iMeynet et all (jl994l ) for masses M > 12 M^, at metal li cities Z = 0.020 and 
Z = 0.040. Note that from H II region measurements, iBresolin fc Kennicuttl (|2002l ) found that the 
metallicity of M83's nuclear region is Z ~ 0.030, which is in between the above two metallicities. 
4) An empirical stellar library of Galactic O and B stars, as well as a new UV high resolution 
theoretical library for massive stars of metallicities Z = 0.013, and Z = 0.030. The difference in 
metallicity between the stellar evolutionary tracks and the theoretical stellar library is due to the 
use of different solar compositions. Indeed, the WM-Basic models use the revised solar composition 
of Asplund et al. (2005), which has a lesser content in metals than the solar composition of the 
stellar evolution tracks. Note that there are no empirical UV stellar libraries for Z > 0.020. 

We describe the empirical and the new fully theoretical stellar spectral libraries in § 13.11 
We discuss the advantages and disadvantages of using the theoretical library in § 13.21 Hereafter, 
we refer to single stellar population (SSP) spectra generated with the empirical stellar library as 
the semi-empirical models, and to SSP spectra computed with the theoretical stellar library, as 
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the theoretical models. We compare theoretical models at Z = 0.020 and Z = 0.040 in § 13.31 We 
compare theoretical models corresponding to different high mass limits of the IMF in § 13.41 Finally, 
we compare the two sets of models in § 13.51 



3.1. Stellar Spectral Libraries 



The semi-empirical models were built from spectral libraries of Galactic O and B stars ob- 
served in the 1000 — 1200 A range at 0.12 A resolution with the Far Ultraviolet Spectroscopic 
Explorer ( Pellerin et al.l I2OO20 . a nd in the 1200 — 1800 A range at . 75 A resolution with the In- 



ternational Ultraviolet Explorer ([Robert et al 



1993 



de Mello et alj|2000i ). Comparisons of the 



semi-empiri cal models wit h observations of local and distant galaxies generally yield quite good 
agreement (|Leithererll2010l ). 



The theoretical models were built from a spectral library of hot massive stars (T j > 20, 000 K, 
M > 8M & ) that was computed as desc ribed in Leitherer et al. (2010), using the stellar atmo sphere 
code WM-Basic ( Pauldrach et al. 200ll ). following the approach pioneered by Rix et al. ( 2004 ). The 
models span the spectral range 900 — 3000 A at high (0.5 A) resolution. This resolution is sufficient 
to resolve the main UV stellar wind lines, which have widths in excess of 1000 km s -1 . 



3.2. New Theoretical Spectral Library 



The strengths and weaknesses of the theoretical stellar library are the following. 1) The two 
main publicly available codes optim ized for modeling the UV spectra of OB stars with winds 
are CMFGEN (jrlillier & Millerl ll998T ) and WM-Basic. WM-Basic, which we utilized, solves the 
radiative transfer and the hydrodynamics self-consistently to derive the density structure of the 
stellar wind, while CMFGEN assumes an analytical approximation for the wind structure a priori. 
2) Wolf-Rayet (W-R) stars, whose stellar evolution is quite uncertain, are modeled using WM- 
Basic atmospheres. This ensures that their continuum fluxes are accounted for, but W-R-specific 
lines are not reproduced correctly. Fortunately, W-R stars in the local universe do not contribute 
significantly to the UV spectrum of OB stellar populations, except for He II A1640, which should be 
viewed with care. In individual W- R stars, the equivalen t width (EW) of the He II 1640 emission 
component ranges from 3 to 100 Al Conti Sz Morris! Il990i ). with a preference for smaller values (3- 
30 Aor so) for late-type WN stars, which tend to dominate observed statistics due to their high 
abso lute visual magnitud e. In starburst clusters, the typical EW of the 1640 bump is about 2- 
3 A (jChandar et al.ll2004l ). If EW=15 A is the typical value for late-type WN stars, this suggests a 
continuum diluting about 80% of the light at 1640 A, and W-R stars cannot contribute more than 
20% to the UV continuum. The dilution is caused by the strong continuum of O main sequence 
and O supergiant stars. 3) We omit the geometrical effect of clumps in the stellar wind, but the 
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associated X-ray emission from shocks propagating in the wind is includeofj]. Clumping affects 
the lines from trace ions, e.g., O VI A1035 and O V A1371 (jBouret et al.ll2005l ). In addition, the 
strength of the wind lines depends on the assumed ionization balance. While photo-ionization is 
the primary mechanism determining the ionization state of the wind, X-rays arising from shocked 
material can be important as well (e.g., iPauldrach et al.l 11994 ). 4) We omit Stark broadening of 
the hydrogen lowest Lyman lines, which is computationally expensive. A pure Doppler profile of 
hydrogen is very similar to one wit h Stark broadening included, as long as the upper principal 
quantum of the transition is small (jRepolust et al.ll2005l ). This is the result of the Stark width's 
dependence on the fourth power of the upper principal quantum number. In practice, this means 
Stark broadening becomes negligible in hydrogen lines at higher energies and originating at lower 
densities. Lyman-a, owing to its large optical depth, is formed far out in the wind even at relatively 
low mass-loss rates. Therefore, Stark broadening is not important at ages younger than 10 Myr. 
5) The evolution of massive stars in binary s ystems and of rotatin g massive stars requires further 
study before it can be included in the models (jVazquez et al.ll2007l ). 6) The models can be used for 
interpreting the rest-frame UV spectra of objects with metallicities other than near-sola r or those 
of the Large (Z = 0.007) and Small (Z = 0.002) Magellanic Clouds (|Maeder et al.lll999l ). 



3.3. Theoretical Spectra at Z = 0.020 and Z = 0.040 



Because there are no empirical UV spectral stellar libraries at Z > 0.020, we have to rely on 
fully theoretical star cluster spectra, at Z = 0.040. Here we discuss differences in the theoretical 
spectra of single stellar populations that were computed with stellar evolution tracks corresponding 
to metallicities of Z = 0.020 and Z = 0.040. 

Line-driven wind the ory predicts lower mass-loss rates an d lower wind terminal velocities for O 
stars of lower metallicity (ICastor et al.lll975l : iPuls et al.lll996l ). and therefore, weaker wind features 
in the spectra of these stars. Observ ations of MS O stars in the Magellanic Clouds and the Galaxy 
support the theory. For instance, see lMokiem et al.l (|2007l ). where empirical stellar wind results for 
the Milky Way, t he LMC and the SMC are compared with the OB star theoretical wind models of 
Vink et al.l (|200lh . which are an extension of line driven wind theory. WM-Basic models do a good 
job of reproducing lower wind terminal velocities at lower metallicity iLeitherer et al.l (|2010i ). as can 
be seen in Fig. \5\ which shows a comparison between the Z = 0.020 and the Z = 0.040 models, 
for ages from one to 20 Myr (as labeled on the right), in the wavelength range 1200 — 1700 A, at 
0.75 A resolution. For instance, compare the absorption components of C IV 1550 at Z = 0.020 and 
Z = 0.040. In addition, the photospheric lines are deeper at the higher metallicity, as expected. 



1 WM-Basic assumes a random distribution of shocks in the wind, where the hot shocked gas is collisionally ionized 
and excite d and emits spontaneo usly into and through an ambient cool stellar wind. This basic model was further 
refined bv lFeldmeier et al.l (1 19971 ) who accounted for the post-shock cooling zones of radiative and adiabatic shocks. 
This approach has been implemented in WM-Basic. 
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3.4. Model Spectra with Different IMFs 

Figure [6] shows semi-empirical and theoretical models for SSPs with different high mass limits 
of the stellar IMF. The figures show ages 2, 3, 5, and lOMyr, and focus on C IV 1550 and Si IV 
1400, whose profiles are sensitive to the distributions of O main sequence and O supergiant stars, 
respectively. All models in Fig. [6] have a lower mass limit of 0.1 Mq. The thin black curves use 
a Kroupa IMF from 0.1 to 100 Mq, while the thick grey curves( red in the electronic version) use 
IMF high mass limits (M up ) of 10, 30, and 50 M . 

Lowering M up , i.e., gradually removing the O and the most luminous B stars, has the following 
effects, (i) As M up decreases, the P-Cygni profiles of Si IV 1400 and of C IV 1550 become weaker. 
In particular, the emission component of these profiles is non-existent at M up = 30 Mq for the 
semi-empirical models, and drops in strength more gradually in the theoretical models. Also the 
strong absorption component of these profiles is purely photospheric at M u = 10 Mq, for both 
models, (ii) For the theoretical models at M up = 10 Mq, the absorption lines are dominated by 
photospheric lines from B stars. 



3.5. Semi-Empirical versus Theoretical Spectra 



Here we discuss differences between the semi-empirical and the theoretical models at Z = 0.020. 
Figure [7] shows the theoretical and semi-empirical models, using the same format of Fig. [5) There 
are several noteworthy features. 1) In the region below 1240 A the semi-empirical spectra suffer 
from Lya absorption in the Galactic plane. 2) The P-Cygni profiles of the resonant doublets of N V 
at 1238.8+1242.8 A (hereafter N V 1240) and of C IV at 1548.2+1550.8 A (hereafter, C IV 1550), 
which originate in the stellar wind, decrease in strength as the O stars evolve and disappear. 3) The 
resonant doublet of Si IV at 1393.8+1402.8 A (hereafter, Si IV 1400) is absent in Main Sequence 
(MS) O stars, and develops P-Cygni profiles in giant and supergiant O stars. This is why the P- 
Cygni profile of Si IV 1400 is strong at ages 3—5 Myr, when O supergiants are present. Note that the 
lack of strong N V 1240, C IV 1550, or Si IV 1400 P-Cygni profiles, in the spectrum of a star cluster, 
indicates either that the cluster did not form very massive O stars, or that the cluster is too old to 
show the signatures of such stars. 4) At 3 Myr, the semi-empirical model has somewhat stronger 
N V, Si IV, and C IV lines than the theoretical model. At this particular time step, stars with 
stronger winds make a stronger contribution to the semi-empirical than to the theoretical spectrum. 
This difference serves as a reminder of a crucial assumption made when linking model spectra to 
stellar evolution models. In the case of the theoretical spectra, one simply adopts the effective 
temperature of the atmospheres for matching the positions of the evolutionary tracks. In contrast, 
for the empirical spectra, we obtain the effective temperatures from spectral types via a temperature 
scale. Star burst99 assumes the r elation between spectral type and effective temperature of O stars 
derived bv lMartins et al.l (120051). W e note tha t switching between t he relations of Martins et al. 
and the earlier work of IVacca et al.l ()1996l ) or ISchmidt-Kalerl (|1982l ) modifies the peak of the Si 
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IV emission by almost 10%. Careful comparisons with observations should provide guidance for 
choosing the preferred approach in the synthesis models. In addition, the WM-Basic models use 
the revised solar composition of Asplund et al. (2005), which has a lesser content in metals than 
the solar composition of the stellar evolution tracks. While the revised solar abundances have a 
small effect on the evolutionary tracks of massive stars, they do affect the computed UV spectra 
via changed opacities and wind properties. 5) Many of the narrow absorption lines in the semi- 
empirical models are of interstellar origin. For the first four Myr, the absorption line Si II A1264 
is stronger in the semi-empirical models because of the interstellar contribution. 6) O V A1371, 
which is strong at ages < 2 Myr, i.e., in the hottest stars, is stronger in the theoretical models due 
to the neglect of wind clumping (jBouret et al.ll2005l ). 



4. ANALYSIS 

4.1. Spectroscopic Ages and Masses 

Table [2] allows an assessment of the uncertainties in the reddened flux at 1500 A, Fi^qq, and in 
the power law index of the reddened FUV continuum, /3 ^ s , for the 13 clusters in our final sample. 
For each spectral exposure and each cluster, columns (2)-(12) in the table give (1) the cluster ID; 
(2) and (3) the right ascension and declination measured from the STIS direct image; (4)- (6) the 
median signal to noise (S/N) of the reddened spectrum; (7)- (9) -F1500; and (10)- (12) /3 ^ s . 

The faint clusters 4 and 14 have the lowest S/N values. The rest of the clusters have reasonably 
good S/N (> 10). The .F1500 fluxes from the slitless exposures taken with the 25MAMA and 
F25SFR2 apertures agree within 10%. Therefore, the Lya airglow has little effect on i^isoo- The 
slit blocks the light from nearby sources and lowers -F1500 by 20 — 50%. The value of Z^^g (i) 
increases by 20% when the wavelength range for fitting the continuum is extended (compare the 
two values in column 10); (ii) varies by at most 10%, due to geocoronal Lya contamination (compare 
columns 11 and 12); and (hi) increases by at most 10% when the slit is removed (compare columns 
10 and 11 or 12). The properties of clusters 1 — 3 derived from the slitless spectra are less reliable 
than those derived from the long-slit spectra, and are only shown for comparison. 

Since the two slitless exposures yield similar values in Tab. El we decided to discard the 
exposure spanning the shortest wavelength range, i.e., the one taken with the F25SFR2 aperture. 
For the nine clusters within the slit, i.e., 1, 2, 3, 4, 10, 11, 12, 13, and 14, the 'good' data span the 
wavelength range 1240 — 1690 A, while for the rest of clusters, i.e., 6, 7, 8, and 9, they span the 
range 1240 - 1670 A. 



We adopted the algorithm developed bv lTremonti et al.l (|200ll ) for deriving the age and mass 



of each cluster from the best fit model to the data. The algorithm takes as input the de-reddened 
cluster spectrum and a grid of Starburst99 models spanning ages from 1 to 20 Myr in step sizes of 
0.1 Myr. The goodness of the fit is characterized by x 2 , where x 2 = (°i ~ m i) 2 Wi/ af , and where 
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Oi represents the observed data for the ith pixel, mi the model data, a the error in the observed 
spectrum, and Wi the assigned weight. The weighting scheme eliminates from consideration the 
interstellar lines, and gives the stellar wind lines a weight 1.5 times greater than that of the con- 
tinuum. Prior to fitting, the resolution of the models was degraded to match the approximate 
resolution of the observed spectra. The spectroscopic masses were derived by (i) fitting power laws 
to the continua of the observed and the model spectra of each cluster, where the model corresponds 
to the best fit age and a stellar mass of 1O 6 M0; (ii) computing the ratio of the observed to the 
model power law at each wavelength point; (iii) taking the average of the sum of the ratios obtained 
in (ii); and (iv) multiplying the result by 1O 6 M0. 



The metallicity of M83's starburst is Z ~ 0.030 with a large error bar (see lBresolin fc Kennicutt 



2002 ) . We considered theoretical models with metallicities of Z = 0.020 and Z = 0.040. 



4.2. Photometric Ages and Masses 



Chandar et al.l (|2010l ) derived the age and extinction for each cluster by finding the minimum 
X 2 when comparing observed UBVIHa magnitudes with those predicted by the Chariot & Bruzual 
(2009, private co mmunication) ste llar population synthesis models, assuming a Milky Way-type 
extinction curve (jFit zpatrick 1999), metallicities of Z = 0.017 and Z = 0.034, and a Chabrier IMF 
from 0.1 — 100 M (jChabrierl 120031 ) . Here, the narrow band filter contains both stellar continuum 



and nebular line emission, i.e. no continuum subtraction was performed for the narrow band image, 
and the predictions were modeled as line plus continuum. 

The mass of each cluster was estimated from its extinction-corrected V band luminosity and 
the age-dependent mass-to-light ratio from the stellar population models, assuming a Chabrier 
stellar IMF and a distance of 4.5 Mpc. The ages and masses derived using this technique are 
compared to the spectroscopic values below. 



5. RESULTS 

5.1. Metallicity and Theoretical Spectra 

On average, the theoretical spectra fit the data better for Z = 0.040 than for Z = 0.020 (on 
average, the reduced x 2 is 1-3 times larger at Z = 0.020 than at Z = 0.040, with a standard 
deviation of 0.3). Note that some of the differences between the data and the best-fit models are 
not due to the metallicity choice, but to the quality of the data, which suffers from crowding, as 
well as from worse spectral resolution than when a narrow slit is used. In any case, on average, the 
spectroscopic ages and masses at the two metallicities are identical. Hereafter, we only consider 
spectroscopic results at Z = 0.020, which is the metallicity of the semi-empirical models. 
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5.2. IMF 

Cluster 14 is the youngest (~3 Myr, average of spectroscopic and photometric ages) and the 
least massive (photometric mass of about 4 x 10 3 M ) in our sample. According to lWeidner et al 



2010 cluster 14 should have a maximum stellar mass of about 80 M (i.e. below 100 M ). Unfor- 
tunately, the spectrum of this cluster is too noisy to reliably pin down its maximum stellar mass. 
Cluster 1 is the next best cluster in our sample to test the upper end of the stellar IMF, and has 
a mean age of ~4Myr and a photometric mass of ~3xlO 4 M . As shown in Fig. [9J left panel, the 
strong Si IV 1400 and C IV 1550 P-Cygni profiles of cluster 1 cannot be reproduced with models 
having maximum stellar masses of 30 M or less. Models with maximum stellar masses of 50 and 



100 M fit cluster 1 equally well. Note that iBresolin &: Kennicuttl (|2002l ) found optical signatures 



of the presence of WN stars in the region of our clusters 1 — 3. The P-Cygni profile of He II at 
1640 A, which is characteristic of WN stars, is barely detected in the FUV spectra of clusters 1 — 3, 
but given the point that we make in § 13.2^ this is not surprising. The rest of clusters with strong 
P-Cygni profiles of N V, Si IV and C IV are consistent with having formed stars more massive 
than 30 M . The right panel of Fig. [9] shows that the upper mass limit of the IMF is harder to 
constrain for cluster 10, which is without P-Cygni profiles. Clusters with strong absorption profiles 
of Si IV and C IV are consistent with having formed at least some B stars, if we assume that the 
absorptions are not of interstellar origin. 

Figures [10] to [12] show the background and reddening corrected rest-frame spectra of the clus- 
ters in our sample, along with best fit semi-empirical and theoretical models corresponding to a 
metallicity of Z = 0.020 and a Kroupa IMF from 0.1 — 100 M . We used the latter models to de- 
termine the spectroscopic ages and masses of all clusters, including those without strong P-Cygni 
profiles of N V, Si IV, and C IV in their spectra. Therefore, we are assuming that clusters without 
strong P-Cygni profiles formed massive O stars but are older than ~12 Myr. Another possibility 
for the latter clusters is that they are younger than ~12 Myr and did not form very massive O 
stars. 



5.3. Spectroscopic versus Photometric Age 

Here, we compare our ages derived using spectroscopic and photometric techniques. We also 
compare our ages with estimates from the literature. 

Our spectroscopic and photometric ages are compiled in columns (l)-(9) of Tab. EJ which 
give (1) the cluster ID; (2) and (4), the spectroscopic ages from the semi-empirical and theoretical 
models; (3) and (5), the reduced x 2 value of the best fits to the observed spectra; (6) and (7), the 
photometric ages from the Z = 0.017 and Z = 0.034 models; (8)-(13) ratios of the age estimates 
derived from the different techniques. The last four rows of Tab. Ogive the minimum, maximum, 
mean, and standard deviation values for relevant columns. Figure [8] shows the photometric age 
determinations for cluster 1 (mean photometric age ~4 Myr) and cluster 10 (mean photometric age 
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~19Myr). They were derived using the technique described in IChandar et al.l (|2010l ). Figure [8] 
also shows the measured V-B and V-I colors for clusters 1 and 10. Considering the mean of the 
spectroscopic and photometric ages, clusters 1 and 10 are representative of young age (< lOMyr) 
and older age (> lOMyr) clusters, respectively. 

The ages estimated from the four different techniques are in excellent agreement. Ages derived 
by comparing UV spectroscopy with semi-empirical and theoretical predictions are within a factor 
of 1.5. This is similar to the la uncertain ties of ~1 Myr fou nd previously for young clusters from 



UV spectroscopy (jTremonti et all l200ll ; IChandar et al.l 120031 ) . Ages determined by comparing the 
HST /WFC3 optical colors of the clusters with predictions from stellar population synthesis models 
at the two metallicities agree within a factor of 2. Finally, the spectroscopic and photometric ages 
are in excellent agreement (within a factor of 1.2 on average). Table [3] shows that most clusters 
are quite young, with ages < 5 Myr. Clusters 6, 7, and 10, which have the oldest spectroscopic 
ages, also have the oldest photometric ages. The worst agreement between the spectroscopic and 
photometric ages is for cluster 8 (difference of a factor of 4). The FUV spectrum of this cluster 
lacks the strong P-Cygni profiles found in young star clusters, yet the cluster has very blue colors, 
resulting in a young (~4 Myr) photometric age. 



Comparing with previous work, our ages agree with the photometric ages derived by lHarris et al 



(|200ll ) based on HST /WFPC2 photometry, except for clusters 6, 7, and 10, which are older than 
6 Myr in our case. We found 3.6, 3.4, and 4.1 Myr, for the ages of clusters 1, 2, and 3 respectively 
(m eans of spectroscopic and pho tometric ages) , which is in excellent agreement with the ages found 
bv iBresolin fc Kennicuttl hood) for region A (all ^4 Myr) using the same spectroscopic data as 
here. On the other hand. IBresolin &: Kennicuttl (|2002l ) determined an age of 7 Myr for their region 
B, from the equivalent width of H/3. Region B corresponds to the location of our clusters 7-9. We 
found mean ages of 15.5 Myr for cluster 7 and 3.4 Myr for cluster 9. As explained in the previous 
paragraph, we were unable to pin down the age of cluster 8, which is also in region B. Finally, 
region D of IBresolin fc Kennicuttl (|2002l ) corresponds to a cluster that we did not include in our 
sample, due to its truncated spectral trace. 

We conclude that our spectroscopic and photometric age estimates are within the expected 
uncertainties for 12 out of 13 of the clusters. The clusters are ~3-20 Myr old and were not all 
formed at the same time. In addition, we do not find any compelling evidence for an age gradient 
along M83's arc-shaped starburst, but rather see clumping of clusters of similar age with very young 
clusters dominating the northern-most and southern-most portion of the starb urst that we s t udied , 
and somewhat older c lusters found between. This is in disagreement with iPuxlev et al.l (|1997l ) 
and iDfaz et all (|2006bl ). who have suggested that there is an age gradient along the starburst, as 
described in the introduction. 
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5.4. Spectroscopic versus Photometric Mass 

Here, we compare the masses estimated for our clusters using spectroscopic and photometric 
techniques. These are compiled in columns (l)-(7) of Tab. [U which give (1) the cluster ID; (2) 
and (3), the masses derived by fitting the semi-empirical and the theoretical models to the FUV 
spectroscopy; (4) and (5) the masses derived from the optical photometry assuming metallicities of 
Z = 0.013 and Z = 0.017; and (6)-(ll) ratios of masses obtained from the different techniques. 

Column (6) of Tab. U] shows that the spectroscopic masses from the semi-empirical models are 
in very good agreement with those from the theoretical models, and that the semi-empirical masses 
are systematically larger by a factor of 1.5 on average than the theoretical masses. The latter is 
because the semi-empirical models have systematically larger values of the mean ratio described in 
(ii) of the last paragraph of § 14.11 Column (7) of Tab. U] shows that the photometric masses are 
not significantly affected by metallicity. Finally, columns (8) to (11) show that the spectroscopic 
masses are systematically larger than the photometric masses by a factor of 4-6 on average. We 
suggest that this is the result of the high background contamination of our wide slit and slitless 
crowded spectra. 

Optical photometry provides more leverage for determining the cluster mass. This is because 
photometry captures light from low and intermediate mass stars, which contribute more by mass 
to the IMF than massive stars. In addition, the total light from the cluster is estimated better 
from the photometry. Therefore, we consider our photometric masses (M p ) as more reliable than 
our spectroscopic masses. 

Our most massive cluster has M p between 2 and 3 x 10 5 M , comparable to t he virial mass 
of the ionizing cluster of 30 Do radus, NGC 2 070 (4.5 x 10 5 M , 1boscIi et al J 12009 ). in the Large 



Magellanic Cloud. According to lLarsenl (|2010l ). young star clusters with masses larger than 10 5 M G 



can last an age comparable to or exceeding the age of the universe. Therefore, this massive cluster 
in M83 could be a globular cluster progenitor, while two other clusters, which have M p ~ 10 5 Mq, 
may also survive. Cluster 14 has M p ~ 4 x 10 3 M and appears to be the most compact in size in 
the FUV image. Finally, the rest of the clusters have M p of a few xlO 4 M . 

Note that our older clusters tend to be more massive than the younger ones. This is probably 
an observational effect. Older clusters at lower masses are harder to detect in the FUV. 



5.5. Intrinsic Reddening 

Columns (l)-(5) of Tab. Ogive (1) the cluster's ID, (2) the power-law index of the continuum 
over the wavelength range 1240 — 1670 A, corrected for a foreground Milky Way reddening of 
E(B — V) = 0.06; (3) the intrinsic reddening of the cluster derived from the FUV spectra; (4) the 
intrinsic reddening derived from the optical photometry; and (5) £i5oo> the intrinsic luminosity of 
the cluster at 1500 A, derived by adopting a distance of 4.5 Mpc to M83. 
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Overall, the agreement between the reddening values derived from the spectroscopic and the 
photometric analysis is fair. For clusters 3 and 14, which have dusty local environments (see Fig. [I]), 
the agreement is good. Indeed, both clusters have large intrinsic spectroscopic and photometric 
reddenings. The agreement is also good for clusters 11 and 12, which appear to be the least 
reddened by both measures of the intrinsic extinction. However, for cluster 1, E(B-V)=0.18 from 
the spectroscopy, while E(B-V)<0.06 from the photometry. Cluster 3 has the largest intrinsic value 
of -L150O) while clusters 4, 13, and 14, which appear faint in the optical images, have the lowest 
values of Li^qq. 



5.6. Diffuse Stellar Field 



In addition to compact star clusters, there is UV emission from the diffuse field star population. 
We obtained the spectrum representing the diffuse stellar field by summing spectra extracted from 
the long slit exposure, at the locations indicated by the dashed lines in Fig. [3l We performed 
no background or reddening correction for the field spectrum. In addition, we did not shift the 
spectral traces in the dispersion direction, due to the lack of narrow interstellar lines that could be 
used for reference. Figure [T3l shows the result. The spectrum lacks the P-Cygni profiles of N V, Si 
IV, and C IV, which are signatures of O stars, but it shows strong absorption components in these 
features, which indicates the presence of B stars. This is consistent with a picture where (massive) 
O stars form mostly in clusters but not in the field, and wh ere the field is dominated by B stars 
from clusters which disru pt rapidly, on timescales of ~10Myr ([Tremonti et al.ll200ll ; IChandar et al 
20061 : IPellerin et ailboOTT ). 



6. SUMMARY AND CONCLUSIONS 



We estimated the ages, reddenings, and masses of 13 clusters within the circumnuclear starburst 
of M83, the nearest such example, based on iKT/STIS/FUV-MAMA imaging and spectroscopy, 
and #ST/WFC3/UVIS photometry. We derived four ages f or each cluster. Th e first two ages 
were obtained by fitt ing the UV spectra with semi-empirical (IRobert et al.lll993l ) and theoretical 



(Leitherer et al 



201CJ) models corresponding to a metallicity of Z = 0.020, giving the most weight 



to the spectral features N V 1240, Si IV 1400, and C IV 1500, which are expected to be prominent 
in the observed spectra of clusters with ages of less than 20 Myr. Note that theoretical models 
with Z = 0.040 yield, on average, ages and masses identical to those at Z = 0.020 . We compared 
our spectroscopic ages with the photometric ages derived by IChandar et al.l (|2010l ) using Chariot 
and Bruzual stellar population synthesis models corresponding to metallicities of Z = 0.017 and 
Z = 0.034. 



The spectroscopic ages derived with the semi-empirical and the theoretical predictions are 
within a factor of 1.2 on average. The spectroscopic and photometric ages agree at a similar level. 



-18- 



Our ages agree with those derived from HST /WFPC2 photometry bv lHarris et aJj (|2Q0ll ) . except 
for clusters 6, 7, and 10, which are older than 6Myr in our case. Our ages for clusters 1- 3 agre e 
with the ages of region A derived from STIS FUV spectroscopy by iBresolin &: Kennicuttl ((2002). 
The clusters are ~3-20 Myr old and were not all f ormed at the same time. We found no ag e 
gradient along M83's starburst, in disagreement with lPuxlev et al.l (|l997l ) and lDiaz et al.l (|2006bl ). 
Our clusters with strong P-Cygni profiles have photometric masses of at least a fewxlO 4 M0, 
seem to have formed stars more massive than 30 M , and are consistent with a Kroupa IMF from 
0.1 — 100 Mq. The rest of the clusters are consistent with having formed at least some B stars. B 
stars dominate the field population. 
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Table 1. OBSERVATIONS 



Type 


Date 


Exp. time 


PA 


Aperture 


AA 


Ext. 






(min) 


(°) 




(A) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


slitless spectra 


1998-05-02 


42 


-24.9 


F25SFR2 


1275-1700 


sfl 


slitlcss spectra 


1998-05-03 


45 


-24.9 


F25SFR2 


1275-1700 


sfl 


slitless spectra 


1998-05-03 


39 


-24.9 


25MAMA 


1150-1700 


sfl 


direct image 


1999-05-19 


23 


32.9 


F25QTZ 


1475-1700 


x2d 


long slit spectra 


2000-04-23 


42 


-9.5 


25MAMA 


1150-1700 


fit 



Table 2. OBSERVED PROPERTIES FROM FUV IMAGING AND SPECTROSCOPY 



ID 


a(J2000) 


<5(J2000) 




S/N a 






Fl500 b 










(h:m:s) 


(":':") 








(erg 


— 1 — t 

s cm 










(1) 


(2) 


(3) 


(4) d 


(5) c 


(6) f 


(7) d 


(8) c 


(9) f 


(10) d 


(H)° 


(12) f 


1 


13:37:00.49 


-29:51:54.93 


18 


17 


27 


6.0E-15 


8.9E-15 


8.0E-15 


-1.3 ( -1.1 ) 


-0.7 


-0.5 


2 


13:37:00.45 


-29:51:55.40 


16 


14 


24 


4.8E-15 


6.8E-15 


6.4E-15 


-1.5 ( -1.2 ) 


-0.9 


-0.9 


3 


13:37:00.54 


-29:51:55.61 


14 


10 


19 


3.8E-15 


4.4E-15 


4.4E-15 


-0.6 ( -0.4 ) 


0.2 


0.8 


4 


13:37:00.42 


-29:51:57.52 


6 






9.2E-16 






-0.6 ( -0.7 ) 






5 


13:37:00.34 


-29:51:58.31 




















6 


13:37:00.37 


-29:51:59.48 




16 


23 




5.3E-15 


5.4E-15 




-1.4 


-1.5 


7 


13:37:00.42 


-29:51:59.83 




24 


35 




1.1E-14 


1.1E-14 




-1.4 


-1.5 


8 


13:37:00.44 


-29:51:59.99 




24 


36 




1.2E-14 


1.2E-14 




-1.6 


-1.6 


9 


13:37:00.43 


-29:52:00.20 




22 


33 




9.5E-15 


1.0E-14 




-1.6 


-1.7 


10 


13:37:01.00 


-29:52:00.71 


21 


22 


34 


7.7E-15 


1.1E-14 


1.2E-14 


-1.5 ( -1.3 ) 


-1.4 


-1.3 


11 


13:37:01.00 


-29:52:01.45 


34 


35 


50 


2.1E-14 


2.4E-14 


2.6E-14 


-2.2 ( -1.9 ) 


-2.1 


-2.2 


12 


13:37:00.54 


-29:52:01.62 


34 


35 


50 


2.1E-14 


2.5E-14 


2.6E-14 


-2.3 ( -2.0 ) 


-2.1 


-2.2 


13 


13:37:00.56 


-29:52:01.91 


22 


23 


35 


9.0E-15 


1.2E-14 


1.2E-14 


-2.3 ( -2.0 ) 


-2.0 


-2.1 


14 


13:37:00.53 


-29:52:02.62 


7 






1.5E-15 






-0.7 ( -0.7 ) 







a Median S/N of the extracted one-dimensional STIS spectrum over the common wavelength range 1300 — 1670 A. 
b Reddened flux at 1500 A measured from STIS spectrum. 

c Power-law index of reddened continuum (F oc A^) over the wavelength range 1300 - 1670 A, which is shared by all 
spectral exposures. For the long slit spectra we also give the value over the extended range 1240 - 1690 A. No data shown 
if/3 obs >0. 

d From STIS long slit spectrogram taken with 1150 — 1700 A aperture. No data shown for clusters outside of slit. 

°From STIS slitlcss spectrogram taken with 1150 — 1700 A aperture. No data shown for clusters with faint spectral 
traces that overlap. 

f From summed STIS slitless spectrograms taken with 1275 — 1700 A aperture. No data shown for clusters with faint 
spectral traces that overlap. 



Table 3. CLUSTER AGES 



ID Spectroscopic 8, Photometric 1 " Ratios from Different Techniques 

(Myr) (Myr) 

semi c x 2d thco e x 2d z017 f z034 s semi/theo z017/z034 scmi/z017 thco/z017 semi/z034 theo/z034 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 



1 


3.9 


3.3 


2.6 


5.4 


4.2 


3.8 


1.5 


1.1 


0.9 


0.6 


1.0 


0.7 


2 


3.4 


4.2 


2.4 


7.4 


4.2 


3.8 


1.4 


1.1 


0.8 


0.6 


0.9 


0.6 


3 


3.5 


4.8 


3.5 


4.0 


4.0 


5.5 


1.0 


0.7 


0.9 


0.9 


0.6 


0.6 


4 


3.5 


1.5 


2.7 


1.9 


4.2 


5.0 


1.3 


0.8 


0.8 


0.7 


0.7 


0.5 


6 


19.8 


1.9 


16.8 


1.7 


10.0 


24.0 


1.2 


0.4 


2.0 


1.7 


0.8 


0.7 


7 


14.2 


2.3 


14.3 


2.9 


9.5 


24.0 


1.0 


0.4 


1.5 


1.5 


0.6 


0.6 


8 


15.6 


2.4 


14.3 


2.4 


4.0 


3.6 


1.1 


1.1 


3.9 


3.6 


4.3 


3.9 


9 


4.5 


5.3 


3.7 


2.3 


3.3 


2.2 


1.2 


1.5 


1.4 


1.1 


2.1 


1.7 


10 


11.8 


2.9 


14.1 


4.7 


13.2 


24.0 


0.8 


0.5 


0.9 


1.1 


0.5 


0.6 


11 


3.4 


6.2 


2.6 


8.8 


6.3 


4.6 


1.3 


1.4 


0.5 


0.4 


0.7 


0.6 


12 


3.4 


7.1 


2.6 


10.4 


6.3 


4.4 


1.3 


1.4 


0.5 


0.4 


0.8 


0.6 


13 


4.2 


2.6 


2.8 


4.9 


9.1 


5.0 


1.5 


1.8 


0.5 


0.3 


0.8 


0.6 


14 


2.2 


2.4 


2.6 


4.1 


3.3 


2.2 


0.8 


1.5 


0.7 


0.8 


1.0 


1.2 


Min 


2.2 


1.5 


2.4 


1.7 


3.3 


2.2 


0.8 


0.4 


0.5 


0.3 


0.5 


0.5 


Max 


19.8 


7.1 


16.8 


10.4 


13.2 


24.0 


1.5 


1.8 


3.9 


3.6 


4.3 


3.9 


Mean 


7.2 


3.6 


6.5 


4.7 


6.3 


8.6 


1.2 


1.1 


1.2 


1.0 


1.1 


1.0 


Stdcv 


5.9 


1.8 


5.8 


2.7 


3.2 


8.8 


0.2 


0.5 


0.9 


0.9 


1.0 


0.9 



a From STIS long slit (25" X 2") 1205.5 - 1690 A spectroscopy for clusters 1-4 and 10 - 14. From STIS slitless 1205.5 - 1670 A spec- 
troscopy for the rest of clusters. 

b From WFC3 optical photometry. 

c Using Z = 0.020 semi-empirical spectra. 

d Reduced x 2 value of the best model fit to the cluster spectrum. 
e Using Z = 0.020 theoretical spectra. 
f Using Z = 0.017 magnitude predictions. 
g Using Z = 0.034 magnitude predictions. 
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Table 4. CLUSTER MASSES 



ID Spectroscopic 3, Photometric 1 " Ratios from Different Techniques 

(M ) (M ) 

semi c theo d z017 e z034 f semi/theo z017/z034 semi/z017 theo/zOI7 semi/z034 theo/z034 

(I) (2) (3) (4) (5) (6) (7) (8) (9) (10) (If) 



1 


1.1E+05 


6.9E+04 


2.5E+04 


2.7E+04 


1.5 


0.9 


4.3 


2.8 


4.0 


2.6 


2 


6.0E+04 


5.2E+04 


1.9E+04 


2.1E+04 


1.2 


0.9 


3.1 


2.7 


2.9 


2.5 


3 


1.4E+05 


1.2E+05 


4.4E+04 


4.7E+04 


1.2 


0.9 


3.2 


2.8 


3.0 


2.6 


4 


3.1E+04 


2.3E+04 


1.7E+04 


1.8E+04 


1.3 


0.9 


1.9 


1.4 


1.7 


1.3 


6 


7.4E+05 


3.8E+05 


6.4E+04 


9.9E+04 


2.0 


0.6 


11.6 


5.9 


7.5 


3.8 


7 


1.0E+06 


6.5E+05 


1.9E+05 


3.1E+05 


1.6 


0.6 


5.5 


3.4 


3.4 


2.1 


8 


9.3E+05 


5.2E+05 


4.1E+04 


3.6E+04 


1.8 


1.1 


22.8 


12.7 


25.6 


14.2 


9 


1.2E+05 


7.7E+04 


1.8E+04 


2.0E+04 


1.5 


0.9 


6.4 


4.2 


6.0 


3.9 


10 


5.7E+05 


4.5E+05 


6.8E+04 


1.3E+05 


1.3 


0.5 


8.3 


6.6 


4.4 


3.5 


11 


9.1E+04 


7.3E+04 


5.0E+04 


5.8E+04 


1.3 


0.9 


1.8 


1.4 


1.6 


1.2 


12 


7.9E+04 


6.3E+04 


5.5E+04 


4.8E+04 


1.3 


1.1 


1.4 


1.1 


1.6 


1.3 


13 


4.3E+04 


2.3E+04 


2.6E+04 


1.6E+04 


1.8 


1.6 


1.7 


0.9 


2.7 


1.4 


14 


3.0E+04 


2.8E+04 


3.9E+03 


4.2E+03 


1.1 


0.9 


7.6 


7.2 


7.1 


6.7 


Min 


3.1E+04 


2.3E+04 


1.7E+04 


1.6E+04 


1.2 


0.5 


1.4 


0.9 


1.6 


1.2 


Max 


1.0E+06 


6.5E+05 


1.9E+05 


3.1E+05 


2.0 


1.6 


22.8 


12.7 


25.6 


14.2 


Mean 


3.3E+05 


2.1E+05 


5.1E+04 


6.9E+04 


1.5 


0.9 


6.1 


4.1 


5.5 


3.6 


Stdev 


3.8E+05 


2.2E+05 


4.7E+04 


8.3E+04 


0.3 


0.3 


5.9 


3.3 


6.4 


3.5 



a From STIS long slit (25" x2") 1205.5-1690 A spectroscopyfor clusters 1-4 and 10-14. From STIS slitless 1205.5-1670 A spectroscopy 
for the rest of clusters. 

b From WFC3 optical photometry. 

c Using Z = 0.020 semi-empirical spectra. 

d Using Z = 0.020 theoretical spectra. 

e Using Z = 0.017 magnitude predictions. 

f Using Z = 0.034 magnitude predictions. 
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Table 5. INTRINSIC REDDENINGS 



ID 




E(B 


-V)i 


E(B 


-V)i 


E(B-V)i 


^1500° 








spec b 


z017 c 


z034 d 


(erg s 1 


A-i) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




1 


-1 


.7 


0. 


18 





.02 


0.06 


1.4E4 


-.38 


2 


-1 


.9 


0. 


16 





.06 


0.10 


9.4E-I 


-37 


3 


-0 


.9 


0. 


28 





.28 


0.20 


2.0E-) 


-38 


4 


-1 


.0 


0. 


23 





.24 


0.20 


5.0E4 


-37 


6 


-1 


.8 


0. 


13 





.20 


0.04 


8.8E-1 


-37 


7 


-1 


.7 


0. 


13 





.22 


0.06 


1.9E4 


-38 


8 


-1 


.9 


0. 


10 





.16 


0.16 


1.5E-4 


-38 


9 


-1 


.9 


0. 


10 





.14 


0.18 


1.3E4 


-38 


10 


-1 


.9 


0. 


14 





.02 


0.00 


1.3E-4 


-38 


11 


-2 


.5 


0. 


05 





.00 


0.08 


1.5E-4 


-38 


12 


-2 


.6 


0. 


04 





.00 


0.00 


1.3E-+ 


-38 


13 


-2 


.6 


0. 


03 





.00 


0.02 


4.8E-) 


-37 


14 


-1 


.1 


0. 


23 





.18 


0.22 


5.7E4 


-37 



a Power-law index of continuum (F oc A* 9 ) corrected for a foreground 
Milky Way reddening of E(B-V) = 0.06, over range 1300-1670 A (slitless 
exposure) or 1300 — 1690 A (long slit exposure). 

Sp ectroscopic red d ening derived adopting the starburst obscuration 
law of ICalzetti et al.l 1120001) and assuming an intrinsic cluster slope of 
P = -2.6. 

c Photometric reddening using color predictions for Z = 0.017. 

d Photometric reddening using color predictions for Z = 0.034. 

e Reddening-corrected luminosity at 1500 A derived from the spectrum 
of each cluster by adopting the spectroscopically derived value of E(B — 
V)i and a distance of 4.5 Mpc to M83. 



Fig. 1.— HST/WFC3/UVIS+STIS/FUV-MAMA composite false-color logarithmic-scale image of 
the two circumnuclear rings of gaseous material embracing M83's starburst. Ha is shown in red, the 
V band in green, and the FUV band in blue. The galaxy's optical nucleus is the large yellow-green 
dot in the upper left. The starburst clusters in our sample are labeled. North is at the top, and 
east is to the left. The location of SN1968L is indicated. 
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Fig. 2.— fiST/STIS/FUV-MAMA data set. Panel (a). 25"x25" logarithmic-scale image (1475 - 
1700 A band), with M83's optical center (cross), compass directions, and spatial scale overplotted. 
Panel (b). 2-D 1150 — 1700 A 25"x2" long-slit spectrogram. The vertical band is due to Lya air 
glow filling the slit. Note the horizontal shift of the spectral traces due to the use of a wide slit. 
Panel (c). 2-D 1150 - 1700 A slitless spectro gram. The Lya air glow contaminates much of the 
image. Note the horizontal shift of the spectral traces due to the use of no slit. Panel (d). 2-D 
1275 — 1700 A slitless spectrogram. Note the lack of Lya air glow contamination. 
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Fig. 3. — Enlarged portions of the STIS images arranged to show the spatial correspondence 
between star clusters and spectral traces. Panels (a) and (c). Image rotated by 42.5° and 57.9° 
East of North, respectively. Lines joining the clusters and their spectral traces are overplotted 
in the two panels. The 2" -wide slit (thick vertical lines) and the extraction regions used for the 
background and the field spectra (solid and dashed horizontal lines within slit, respectively) are 
overplotted in panel (a). Only the top background region is shown in panel (c) since the bottom one 
falls outside of the figure. Panels (b) and (d). 2-D 1150 - 1700 A long-slit and 1275 - 1700 A slitless 
spectrograms. The size of the box used for extracting the 1-D spectra is shown. 
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Fig. 4. — Plots showing the spatial correspondence between the star clusters in the STIS image, 
and the long-slit and slitless spectral traces (left and right panels, respectively). In each panel, 
the top curve represents the cross dispersion profiles of the spectral traces, while the bottom curve 
represents the radial profiles of the clusters. The clusters in our sample are labeled 1-14. The 
background subtraction regions are labeled "sky". The dotted curve is the profile of a point-like 
white dwarf observed with the same configuration that was used for the long-slit spectrum. The 
profiles were obtained and scaled as explained in the text. 
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Fig. 5. — Comparison of SSP theoretical spectra with Z = 0.020 (thin black curves) and Z = 0.040 
(thick grey curves, red in electronic version). A Kroupa IMF from 0.1 to 100 Mq was used. Each 
spectrum represents an age step from 1 to 20 Myr, as labeled on the right. The spectral features 
discussed in the text are labeled at the bottom. 
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Fig. 6. — Comparison of SSP spectra with high mass limits of 10, 30, 50, and 100 M & for the 
IMF. The thin black curves correspond to a Kroupa IMF from 0.1 to 100 Mq, while the thick grey 
curves (red in electronic version) correspond to the high mass limits indicated in each panel. Semi- 
empirical and theoretical models are shown on the left and right panels, respectively. Z = 0.020. 
The region 1350 — 1650 A is shown at ages 2, 3, 5, and 1 Myr, as labeled on the right. The age 
sensitive features Si IV 1400 and C IV 1550 are labeled in the bottom two panels. 
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Fig. 7. — Comparison of SSP spectra generated with theoretical (thin black curves) and empirical 
(thick grey curves, red in electronic version) stellar libraries. Z = 0.020. A Kroupa IMF from 0.1 
to 100 Mq was used. Each spectrum represents an age step from 1 to 20 Myr, as labeled on the 
right. The spectral features discussed in the text are labeled at the bottom. 
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Fig. 8. — Photometric measurements and age determinations for clusters 1 (left) and 10 (right). 
The top panels show the predicted SEDs from Chariot & Bruzual models for Z = 0.017 (2009, 
private communication), for the age and extinction that best fit the observations (solid line). The 
data points represent measurements in U, B, V, Ha, and I filters, and have been corrected for 
extinction in both the Milky Way and in M83 (values for the latter are listed in Tab. [5]). Note 
the strong nebular emission in the Ha measurement for cluster 1. The bottom panels compare the 
U — B vs. V — I colors for each cluster with model predictions for metallicities of Z = 0.017 (solid 
line) and Z = 0.034 (dashed line). The arrow shows the direction of reddening. The observed 
cluster colors have been corrected for reddening in the Milky Way but not in M83. The triangles 
mark ages of 1, 10, 100, and 1000 Myr along each model, starting from the upper left. 
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Fig. 9. — Observed spectra of clusters 1 and 10 (thin black curves) versus semi-empirical models 
corresponding to IMF upper mass limits of M up = 10, 30, 50, and 100 M Q (thick grey curves, red 
in electronic version). 
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Fig. 10. — Background and reddening corrected rest-frame spectra of clusters 1 — 4 (thin black 
curves), with best fit models overplotted (thick grey curves, red in electronic version). The left 
panels show the semi-empirical models while the right panels show the theoretical models. The 
y-axis gives the luminosity as a function of wavelength from ^1200-1700 A. 
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Fig. 11. — Same as Fig. [TUlbut for clusters 6 — 9. 
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Fig. 12. — Same as Fig. [10] but for clusters 10 — 14. 



-38- 




Fig. 13. — Extracted spectrum of the diffuse stellar field. For this spectrum, we performed no 
background subtraction, reddening correction, or horizontal shifting of the spectral traces. 



